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Abstract 

We introduce in this paper an elliptic dynamical reflection algebra describing an 
SOS model with reflecting end. Using factorizing Drinfel'd twist, we compute 
the partition function of this model with domain wall boundary conditions. We 
show that it can be represented in the form of a single Izergin determinant. 
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1. Introduction 

SOS type models in statistical mechanics play an important role for math- 
ematical physics. Since the pioneer works of Baxter on the eight vertex model 
[2], they arise as a necessary step towards resolution of vertex model without 
charge conservation [l(| . As they are described by dynamical Yang-Baxter 
algebra, they also arise as underlining algebraic structure of conformal field 
theories [l5j |. From the algebraic point of view, they are intensively studied 
due to their relations to quasi- hopf structures Q and current algebras 20l. 22\. 



Furthermore, they are related to combinatorics and dynamical enumeration of 



alternating sign matrix 2l|. 



Within the Quantum Inverse Scattering frameworks, partition functions are 
essential for the study of correlations functions @, 0, [Til El- The partition 
function for such elliptic SOS model with domain wall boundary conditions 



has already been computed by different methods in [22J, |21|, but due to the 
dynamical nature of the underlining algebras, it is presented as a sum of matrix 
determinants. Recently, partition function of an alternative trigonometric SOS 
model with reflecting end has been computed as a single determinant [l7| . It is 
shown there that such boundaries permit to avoid the difficulty of the inherent 
dynamical algebra, by means of a dynamical reflection algebra. We generalize in 
this letter the result to the elliptic case, which is the most general SOS models, 
by use of the concept of factorizing Drinfel'd twist. 
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2. Preliminaries: elliptic theta functions 

Let r be a fixed complex parameter such that: Im(r) > and denote 
p = e T . Throughout this paper, we will use the notation 



,2A\ 



(1) 



i=0 



Up to a multiplicative factor, h(X) equals the Jacobi theta function 9i(i\) [251 ] - 
This function is odd and satisfy the addition rule 

h(x + u)h(x — u)h(y + v)h(jj — v ) — h(x + v)h(x — v)h(y + u)h(y — u) (2) 
= h(x + y)h(x — y)h(u + v)h(u — v). 

In the degenerate case, we have: lim p _>o h(X) = 2sinh(A). 

Definition 2.1. f is a theta Junction of norm t and order N if there exist N 
constants {£i}i=i,...,Ar cind £i = ^ such that 



N 



f(\) = l[h{x + i l ). 



(3) 



i=i 



Using this definition, we have the classical theorem 24 1 
Theorem 2.1. Let f be 

N 

j i=l 



(4) 



where for any j , y\ —1 £f — t. Then f is a theta function of order N and norm 
t. 

This theorem gives a powerful tool to prove theta functions identities. 



3. Algebraic framework 

We start by introducing the dynamical reflection algebra which is underlin- 
ing integrability of the elliptic SOS model with reflecting end as we will show in 
this paper. This algebra is actually built as a comodule over the elliptic quan- 
tum grou p E misl^]- The main object for defining the elliptic quantum group 
E T , v (sl 2 ) lH is the dynamical R matrix, R : C x C — > End(F <g> V), V - C 2 



R(X;9) = 



\ 



(RXX(X;6) 

R+Z(X;9) R++(X;9) 

R+±(X;0) RZ+(X;0) 

V RZZ{X;9)J 



(5) 
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which satisfies the dynamical Yang-Baxter equation 



iiia(Ai - \2\0- w z 3 )R 13 (\ 1 - X 3 ;9)R 23 (X 2 - X 3 ;9- w \) 

=R 23 (X 2 - \ 3 ;0)R 13 (\ 1 - A3; - t^)J2ia(Ai - A 2 ; 0), 

(6) 

where we denote by a^' v ' z the usual Pauli matrices in the two dimensional s pac e 
V a ~ C 2 . We are interested here in the elliptic solution of this equation |19j . 
which is the most general case 



Rli(\;8) = RZZ{X;9) = h(X + r)) 

h(X)h{ 
K 

h(r])h(6 - A) 



RtZ{X;6)=RZi(X;-B)= h ^-^ (7) 



R±+{\;0)=R+±(\;-e) 



HO) 



The elliptic E T ^{sl 2 ) quantum group is the algebra generated by meromorphic 
funtions of the generator of f), the Cartan subalgebra of sl 2 , that we denote 

by a z , and the matrix elements of T(A; 6) = ( 2 nnlel I G End(C 2 ) 



\C(X;6) D{\;0) / 

with non-commutative entries, satisfying the dynamical Yang-Baxter algebra 
relations 



Ri2{Xi -X 2 ;9- W Z )T 1 (X 1 ;0)T 2 (X 2 ;6 - W {) 

=T 2 {X 2 -9)T 1 {X 1 -9- w z 2 )R 12 {X 1 -X 2 -e). (8) 

We are interested here only in diagonalizable ()-module V where the weight 
zero property holds 

[T (X;9),a z +a z v ] = Q. (9) 

We define the following Dynamical Reflection Algebra generated by mero- 
morphic functions of a z E f) and the matrix element of T(A; 9) £ End(C 2 ) with 
non commutative entries subject to the relations 



Ru(Xi - A 2 ; - r 1 a z )Ti(X 1 - 1 9)R 21 (X 1 + X 2 ;9- w z )T 2 (X 2 ; 6) (10) 
= r 2 (A 2 ; 9)R 12 (X 1 + X 2 ;9- W Z )T 1 (X 1 ;9)R 21 (X 1 -X 2 ;6- rya 2 ). 

Let /C : C x C — > End(C 2 ), be a (scalar) representation of this algebra in 
C (i.e C-number matrix) viewed as one dimensional f)-module of sl 2 with the 
standard action on v G C, a z .v = 



3 



Rn(Xi - A 2 ; 0)£i(Ai; 9)R 21 (X 1 + A 2 ; 0)/C 2 (A 2 ; 0) 

=/C 2 (A 2 ; 6)R 12 (\i + X 2 ;9)tC 1 {X 1 ;9)R 2 i{X 1 ~ A 2 ; 6). (11) 

This is essentially the reflection equation introduced in with the dy- 
namical i?-matrix instead of the usual one. A representation as above is said of 
weight zero if 

{JC Q (X;9),a z }. (12) 

This means that JC is diagonal solution of the above equation. Let T(X; 9) a 
weight zero representation of the dynamical Yang-Baxter algebra in V . Then 

T(A; 9) = T(A, 9)JC(X; 0)T- x (-A; 9), (13) 

is a weight zero representation of the dynamical reflection algebra in C ® V. 

In this letter, we are mostly interested in a particular representation of this 
dynamical reflection algebra, which is built on the well-known evaluation rep- 
resentation of E T:Tj (sl 2 ) in the space V = (Si^Cf which is constructed from the 
i?-matrix jj]) 



N 

T (\; 9) = RmiX-h; d-rlY,*!) -Ron (X-Zn;0) (14) 

i=2 

_(A(X;9) B(X;9y 
\C(X;6) D{X-9) / 

and the specific diagonal solution of (fTTT) 

/ h(8+(-\) n \ 

me)= TOAJ , (15) 
\ u MC+A)/ 
which depends on an arbitrary parameter £. So our main object of study is the 
boundary monodromy matrix, a representation of (fTOj) in C x V, V = ®£iC? 



T(A;0) =7(A)T(A;^)X;(A;e)T- 1 (-A;^ (16) 

M(A; 0) B(A;0)\ 
~^C(A;0) 2)(A; 0) J ' 

with normalization coefficients 

N 

7(A) = (-1) W J] h ( X + d - VMX + & + rf). (17) 

i=i 

This representation has a clear statistical mechanics interpretation, it de- 
scribes an elliptic SOS model with reflecting end. Before introducing this sta- 
tistical mechanics model, we give in the next section convenient expression for 
the boundary monodromy operators (|16[) . 
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4. Symmetry and Drinfel'd twist 

In this section, we collect some useful relations of the i?-matrix and corre- 
sponding monodromy operators. This relations will help us to obtain manage- 
able expressions of the boundary operators (IT6l) with the help of the factorizing 
Drinfel'd twist. 



4-1. R-matrix and boundary operators 

The i?-matrix satisfies three important properties 

1. Zero weighfl 

[at+al,R 12 (X;e)} =0 (18) 

This symmetry reflects the six vertex texture of the statistical weights: 
R 1 ^ = unless a + (3 = fi + v. It is easy to see that this relation induces 
a similar relations for the transposed R matrix 

[af-af,J&(A;0)]=O (19) 

2. Unitarity: 



#i 2 (A;0)-R 2 i(-A;0) = -h{\ - rj)h(\ + ij)Id (20) 
3. Crossing Symmetry: 

We write the crossing relation for the dynamical i?-matrix in the following 
compact form: 

- a\ : R%{-\ - m 6 + V al) ■ °l = ^i(A; 9) (21) 

where we assume the following normal ordering : the erf in the argument of the 
i?-matrix (which does not commute with it) is always on the right of all other 
operators involved in the definition of R. 

4-2. Boundary monodromy matrix 

Using the crossing relation (f2~Tj) and the zero weight for the transposed R- 
matrix (fT9|) we get 

JV 

7(A)T-VA; 9) =R m (X + £ N ; 9)...R 10 (X + 6; - r?X>i ) 

i=2 

= 7 (AK¥»(-A- I) ;» + ^K 
h(9 -t]S z ) 

h(9) ' (22) 



1 This property is sometimes referred in statistical mechanics as the Ice Rule 
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with normalization coefficients 

7 (A) = (-If, (23) 

and the short notation: Oy = S z = u t ■ A very important decomposition 

of the B operators is given by means of the generalized crossing relation for the 
bulk monodromy matrix (|22[) . It implies for the B operators 

B(X; 9) = 7 (A) (fCZB(X; 9)A(-X - 77; 9 + rj) - K%A(X; 0)B(-\ - 77; 6> - rfy 
h(9-nS z ) 



h{6) 



(24) 



Furthermore, the equation (|10p contains the commutation relations for the 
generators: A(X;9), B(X;9) : C(X;9) and T>(\;9). The only one which is im- 
portant for the computation of the partition function is the relation for the B 
operators 

B(\ l ;9)B{\ 2 - 9) = B{\ 2 - 9)B(X 1 ;9) (25) 

4-3. Symmetric representation and the F-basis 

Representation of Drinfel'd twist Q was first applied by Maillet and Sanchez 
de Santos [l3| in order to obtain completely symmetric representation of the 
bulk monodromy operators for Yang-Baxter type algebra, which are highly non 
local in terms of the quantum local operators. The idea is to perform a change 
of basis in the space of states where the bulk monodromy operators remains 
completely symmetric. 

This representation is very useful for the computation of partition functions, 
scalar products and correlation functions of integrable spin chains with periodic 
boundary conditions [8|, |9| but also for diagonal boundary conditions |lO|,lll|, as 
it reduces drastically the combinatorial difficulty of handling highly non-local 
representation. 

In this representation, not only the bulk monodromy operators remain sym- 
metric, we also have a direct insight into the analytical property of this opera- 
tors. 

For the evaluation representation of dynamical Yang-Baxter algebras, such 
representation of the bulk monodromy operators was proposed in [lj]. This 
construction is based on a dynamical _F-matrix which factorizes the dynamical 
i?-matrix in the following way 

F 21 (-X;9)R 12 (X;9) =F 12 (X;9). (26) 

After a suitable co-product over all quantum spaces it leads to a change of basis 
Fuy where the bulk operators A(X; 9), B(X; 9) (fTl"]) have symmetric expressions 

A(X; 9) - F {£} (6)A(X; 0)F^\ (9 - r,) (27) 

h(9- v ) ^ N {h(X-ti + ri) 

h(9 + r)( 
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and 



B(\;6)=F {6} (6)B(\;6)F-](8 + r,) 



(28) 







A very important property is that the reference states |0) = n*Li t& an d 



|0) = YiiLi l-Si are ^ e ft an d right invariant under the action of 

2W0)|O> - FZ\(9)\0) = \0)A0\F {(} (9) = (0\FzU6) = <0|, 



(29) 



and 



F {6} (8)\0) = F { "g(0)|O> = <0|, (O\Fm(0) = (0\F£(6) = <0|. (30) 

Using the decomposition (|24|) of the boundary operator B(X; 8), it is easy to 
compute its expression in this new basis 



B(X;9) = F {6} (9)B(X;9)F^(9) 



(31) 



N 

KA)E{ 



i=l 



h(9 + { + X) h(( + x) 



h(2X)h(r}) 



fh(X + ^)h(X-^+T,) 

>iV I n h(A-£,- )ft(A+gj +rj)fe(g, -gj +r,) 



fe(6>-7?S z ) 



5. Elliptic SOS model with reflecting end 

Let us now introduce the SOS model, which is a two-dimensional statistical 
mechanics lattice model where Boltzmann weights are attached to each face. 
There are 6 possible face configurations where each height 9 can differ only by 
±i] for adjacent sides 





9-27] 


8 + 1] 


8 + 2 v 


8-i] 


8 


9 


8-1] 8 


8 + t] 8 


8 + T] 



8 + t] 


8 8 + i] 


8 8-i] 


8 


8 


8-i] 8 


8 + i] 8 


8-i] 
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/C+(A;6>) fCZ(\;6) 

Figure 1: Boundary configuration with external height 8. 



The corresponding statistical weights R^ are collected into the dynamical 
R matrix (JTJ). In the limit p — > 0, the model becomes equivalent to the trigono- 
metric SOS model. Furthermore, the limit 9 — > oo gives rise to the well-known 
six vertex model |12j . Dealing with boundary weights requires to introduce the 
boundary matrix K.(X; 8) which we choose to be (fT5|). Notice that in the trigono- 
metric case, /C(A; 6) reduces to the trigonometric solution of [l7j . We consider 
this model with a reflecting end, which means that each horizontal line makes 
a U-turn on the left side of the lattice. As we choose a diagonal solution, it 
produces two configurations characterized by the weig hts AC|(A;0). 

It is important to note that such reflecting end imposes a constant external 
height 9 for the left side of the lattice. This reflecting end leads to different 
paramctrizations of the weight if they are in the two different half the rows. 
Indeed, parametrization should respect row and line multiplication, and also 
some fundamental symmetry of the R matrix (as the zero weight (|18p V We can 
easily check that this convention leads to a well defined inhomogeneous model: 



T>c— cZ,a — c/ \ 
H a-b,b-d\ A 3 



a 


b 


c 


d 


e 


f 



R 



f-d,e-f 



(A,--6;e) 



Here we consider the model with domain wall boundary conditions, the 
heights decrease from left to right on the upper boundary, the heights grow 
from left to right on the lower boundary. As left external height is fixed these 
two conditions determine completely the configuration on the right boundary 
(heights decreasing in the upward direction). 
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i 




9-(N-l)r| 


8-Nt] 








9-(N-1)ti 








9-(N-2)t) 








9-(N-3)ti 






















9+(N-3)t] 








9+(N-2)r| 








0+(N-l)T| 


e+ri 




9+(N-l)r| 





Figure 2: Domain Wall Boundary Conditions 



6. Partition Function 

The partition function of the SOS model introduced in the previous section 
can be written in terms of the boundary monodromy matrix 



N 



N 



N 



N 



N 



A' 



l[B(K-9)\0) 

i=i 

N 



(32) 



We follow the standard way to compute the partition function [12|, |6j , first we 
establish a set of properties defining it in an unique way and then we propose a 
determinant formula which satisfies all these conditions. 

The partition function (|32|) satisfies the following properties: 

i) -Z/v,2Ar({A}, {£}, 9) is symmetric in X{. This property follows from the com- 
mutation relation (|25|) for the operators B. 



ii) Zm,2n{{^}^ {£,}, 9) is symmetric in £j. This is a direct consequence of the 
Dynamical Yang-Baxter Equation ([5]). 

It is sufficient to insert — — ^^+2) m f32|) to get the 

symmetry for any elementary permutation £j O 
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iii) For each parameter A^ the normalized partition function 



N 



z N MW, {*}, &) = II h{9 + C t { 2x^ + Xi) zN ' 2Ni{xh U} ' (33) 

is a theta function of order 2N — 2 and norm (N — l)r) with respect to 
the variable Xi. To prove it, we use the representation of the partition 
function in the F-basis and consider the action of the most right B(Xn', &) 
operator. Due to the symmetric representation of B(Xn', 0) in the F-basis, 
the operator 

M2A^ B{Xn ' 9) 

acts due to (f3"Tl) as 

N N 

which is a theta function of the desired form thanks to the standard theorem 
of Section 1. This property remains true for any Xi with the help of i). 

iv) For N = 1 the partition function is just a sum of two terms 

h(6 + C - A) 

h(e + <; + x) 

, ^(C - A) 
h(C + A) 

as there are only two configurations possible. 

v) Recursive relations. 

There are two points where we can easily establish recursive relations, fixing 
the configuration in the lower right or the upper right corner by setting 
Ai = £i or Ajy = — £i- It is easy to see that it leads to the following 
recursive relations 



h(x-c)h(e + x + c) 

h{X + £,)h{6 - X + ) (34) 



Zn,2n ({A}, 



= h( V )h(C - Ax 
MC + Ai) 



n,,, /■ \ h(6 + (N - 2i)r)) 

N 

x Y[ h(\i - f< + ??)/i(Ai + & + r?)/i(Ai - 6 + »?) 

X V-l),2(iV-l)({A}2...JV I Uk.JV, {0}) (35) 
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Zn,in ({A},U},W) 



= h( v )h(0 + (-\ N ) 
h{8 + ( + \ N ) 



z=l 

AT 

X 
i=2 

x V-i),2(^-i)({A}i...jv-i, UK.jv, {#}) (36) 



/i(0 + (N - 2i + l)?j) 
] { ^(Aw + & + vMXn ~ & + v)H*i-i + Ci + »?) 



Lemma 6.1. T/ie sei of conditions i)-v) uniquely determines the partition June- 
tion Z N . 2N ({\},{£},{9}). 



Indeed, it is sufficient to observe that the normalized partition function (|33[) is 
a theta function of order 2N — 2 and norm (N — l)rj in each parameter A^. So 
we need 2N — 1 independent conditions to uniquely determine it. Using the 
symmetry ii) the recursion relations v) can be established for any point A; = £j 
( or A; = — £j ). Hence we can prove by induction starting from the case N = 1 
that the partition function is uniquely determined as we need. 

Theorem 6.1. The partition function of the elliptic SOS model with reflecting 
end can be represented in the following form 



ZnMW, {£}, M) = 7(A) det M 13 ]J 



' h(9 + r](N -2i)) 



h(6 + 7](N - 2i + 1)) 



TV 

n + &X*i - Zj)h(\i + & + v)HXi - & + n) 

x !___ 

II KZi + &)KSi - Zi)H*j - A i )/i(A j + K + rf) 

\<i<j<N 

where the N x N matrix Mij can be expressed as 

= h(9 + ( + Q h(C-Q 
hJ h(9 + ( + X t ) h{( + \ t ) 



(37) 



h{2X l )h( v ) 



h(\i - + r))H*i + + vM^ - Zj)h(\i + £,) 
To prove the theorem, it is sufficient to check the property i) to v). 



(38) 



Remark 6.1. The partition function reduces to that of the trigonometric SOS 
model with reflecting end in the limit p — >• 
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Taking the limit of (f3"T| we obtain up to irrelevant numerical factor 

iV 

n sinh(Ai + £j) sinh(A i - £,■) sinh(Aj + £j + 77) sinh(Ai - + 77) 



I~[ sinh(^j + sinh(^ — sinh(Aj — Xi) sinh(Aj + Aj + 77) 

l<i<j<N 

(39) 

where the N x N matrix M^ is 

_ sinh(g + C + 6) sinMC-^j) 
lJ sinh(0 + C + A,-) sinh(C + Ai) 

x sinh(2Ai)sinh?7 

sinh(Aj — + 77) sinh(Aj + £j + 77) sinh(Ai — £,) sinh(Aj + £j) ' 

We recover in this way the previous result of [l7j |. 
Conclusions 

In this paper, we have computed the partition function of the elliptic SOS 
model with domain wall boundary conditions and one reflecting end. The main 
result is that this partition function is expressed as a single determinant, gener- 
alizing the recent observation in the trigonometric case that SOS type models 
with such boundaries give rise to a simpler representation of their partition 
functions. Recently, the trigonometric version of this SOS model was associated 
to open XXZ spin chains with particular boundary conditions in a very natural 
way by means of gauge transformation (Vertex-IRF, [18|). It seems to be natu- 
ral to ask for such relation between the result of this paper and the open XYZ 
spin chains, leading to the possibility to study XXZ and XYZ spin chains with 
boundary in a unified framework in opposition to the periodic case. 
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